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ABSTRACT 

We performed a series of hydrodynamical calculations of an ultra-relativistic jet propagating through 
a massive star and the circumstellar matter to investigate the interaction between the ejecta and the 
circumstellar matter. We succeed in distinguishing two qualitatively different cases in which the ejecta 
are shocked and adiabatically cool. To examine whether the cocoon expanding at subrelativistic speeds 
emits any observable signal, we calculate expected photospheric emission from the cocoon. It is found 
that the emission can explain early thermal X-ray emission recently found in some long gamma-ray 
bursts. The result implies that the difference of the circumstellar environment of long gamma-ray 
bursts can be probed by observing their early thermal X-ray emission. 

Subject headings: gamma-ray burst: general - radiation mechanisms: thermal - shock waves - super- 
novae: general 
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1. INTRODUCTION 

Since the discovery of gamma-ray bursts (GRBs), nu- 
merous studies have been done to understand their pro- 
genitors, the mechanism to produce their hig hly energetic 
emission, and t he central engine (see, e.g., [1999; 
iMesza ros 2006, for review). It is currently known that 
long GRBs are triggered by the gravitational collapse of 
massive stars. The spatial an d temporal coinciden ce of 
GRB 980425 and SN 1998bw (|Galama et al.l Il998) has 
revealed the connection between long GRBs and a spe- 
cial class of type Ic supernovae (broad lined type Ic 
SNe), i.e. , the firmly established SN-GRB connection 
(see, e.g., Wooslev & Bloom 2006). For example, well- 
known GRBs associated wi t h SNe are GR B030329/SN 
2003dh (jHiorth et aLMgOOgHS tanek et al. 200 31), GRB 
0602 1 8/SN 2006aj (ICampana et al.l 120061 : iPian et al.l 
2006: iMazzaliet all 12006ft. GRB 10Q316D/SN 2010bh 
" "ano et al.l 120111 : iBufano et al.l 120121 : lOlivares E. et al.l 
120121 ). The increasing number of detected samples of 
GRB-associated SNe has enabled us to investigate their 
circumstellar environments. 

Especially, whether the circumstellar matter (CSM) of 
the progenitor is dilute or dense is of particular inter- 
est, because it is expected that the CSM interacts with 
the ejecta and results in producing high-energy emission. 
The CSM may originate from the stellar material ejected 
prior to the explosion as a wind or the common envelope 
if the progenitor of the G RB was in a binary system 
(jPodsiadlowski et al.l I2004D . 

Recently, it is reported that thermal components are 
found in X-ray spectra of some long GRBs, which are 
taken by Swift satellite 100-1000 seconds afte r the trigger 
(ICam pana et al. 2006; Starli ng et all 120111: iPage et al.l 
120111 : IStarling et al.ll2012USparre fc Starlingl 120 12ft . The 
component is seen as an excess superposed on a power- 
law non-thermal component that is usually attributed to 
synchrotron emission from the forward shock, i.e., the af- 
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terglow emission. Spectral analyses reveal that the com- 
ponent can be fitted by a single blackbody spectrum with 
temperature of k^T — 0.1-0.9 keV (see, IStarling et al.l 
2012). The luminosity ranges from 10 45 to 10 49 erg s 1 . 
The contribution of the thermal emission to the total 
X-ray flux is typically a few % up to several 10 %. The 
emitting radii inferred from the fitting results are 10 12-13 
cm, which are much larger than the typical radius of 
the progenitor star < 10 11 cm. Their durations are 
several 100 seconds, up to 1000 seconds for the longest 
case, GRB 060218, which is classified as a low luminosity 
GRB associated with a supernova SN 2006aj. The num- 
ber of GRBs whose spectra exhibit the thermal compo- 
nent now reaches severa l dozens (see. IStarling et aO2012t 
ISparre fc Sta7Iinll20T2h . 

Several models to explain this emission component 
have been presented. As an example, it is proposed that 
the supernova shock breakout c an be responsible for the 
emiss ion of some GRBs (e.g., IWaxman et al.l 12 0071 : lO 
12001 . On the other hand, for GRB 060218/SN 2006aj, 
it is pointed out that the radiated energy and the in- 
ferred emitting radius are too large to ascribe the emis- 
sion to the s upernova shock breako ut from the progenitor 
surface (e.g., Ghisellini et al. 2007). Therefore, some au- 
thors ascribe the large emitting radius to the presence of 
a stellar wind with a high mass-loss rate. In this model, 
the shock emerges from the photosphere located in the 
wind. Another proposed model is the cocoon emission. 
The cocoon is a hot plasma resulting from the interaction 
between the jet and the stellar material. It emerges from 
the star at the same time the collimated jet penetrates 
the stellar surface and then expands s pherically at mildly 
relativist ic speeds. iPe'er et all (|2006f ) investigated emis- 
sion from the cocoon by combining a numerical radiative 
transfer calculation with an analytical treatment of the 
dynamical evolution of the cocoon. While their model 
is easy to treat, it is necessary to check whether some 
parameters used there, such as, the total energy of the 
cocoon, are realized in actual situations by using hydro- 
dynamical calculations. In particular, by performing hy- 
drodynamical calculations, one can estimate the amount 
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Fig. 1. — Color-coded Lorentz factor and density distributions at t = 6 (top left), 10 (top right), 100 (bottom left), and 10 3 (bottom 
right) s for the model with M = 10 -7 M© yr _1 . 



of energy deposited into the cocoon out of the total in- 
jected energy in a self-consistent way. Furthermore, the 
large emitting radii inferred from spectral analyses indi- 
cate that the emission comes from the region where the 
CSM is expected to be present. If so, the ejecta-CSM in- 
teraction must give rise to thermal X-ray emission. This 
effect should also be investigated by hydrodynamical cal- 
culations. In addition, the cocoon emission might be im- 
portant as a source of seed photons for inverse Compton 
to produce high-energy photons with ener gies of ~ 100 
MeV, as pointed out bv lToma etall (|2QQ9h . 

In this Letter, to investigate the interaction between 
the ejected matter and the CSM, we perform special rel- 
ativists hydrodynamical calculations of the propagation 
of a relativistic jet emanating from a massive star in the 
CSM. In Section 2, we describe our method to calculate 
the evolution of the jet and the interaction with the CSM. 
Results of the hydrodynamical simulations and the ex- 
pected light curves of the emission from the cocoon are 
presented in Section 3. Finally, in Section 4, we discuss 
implications from the results and conclude this Letter. 

2. METHOD 

In this section, we briefly explain setups of the hy- 
drodynamical calculations performe d in this study . The 
detailed code description is found in iSuzukil (j2012h . 

2.1. Hydrodynamics 

We perform hydrodynamical calculations of the prop- 
agation of an ultra-relativistic jet in a massive star and 
the subsequent interaction with the CSM by using the 
special relativistic hydrodynamics code in 2D spherical 
coordinates (r, 0) developed by one of the authors. In this 
code, we adopt a mapping procedure, in which the width 
of the radial zones is doubled as the jet head reaches a 



fraction (~ 0.9) of the maximum of the radial coordi- 
nate, in order to calculate the propagation of the jet till 
t ~ 1800 s. Thus, the radial resolution becomes coarser 
as the time elapses. At t = 0, the radial coordinate 
ranges from r = 10 9 cm to r = 10 11 cm. At the end of 
the calculations, the maximum of the radial coordinate 
reaches r ~ 6 x 10 13 cm. The radial zone is divided into 
N r uniform cells and the number N r — 1024 is fixed. 
The angular coordinate ranges from 6 = to 6 = tt/2 
and is composed of Nq = 256 uniform cells. 

2.2. Simulation setup 

As a presupernova m odel, we adopt 16TI model in 
IWooslev fc H eger (2006), which is commonly used in cal- 
culations of collapsar jets. In this study, we consider 
several models to clarify the effect of the ejecta-CSM 
interaction. Since the spatial distribution of the CSM 
is highly uncertain, we adopt the simplest steady wind 
model whose density profile is given by, 

The density profile is uniquely determined for a given 
ratio of the mass-loss rate M and the wind velocity 
v w . In this study, the wind velocity v w is fixed to 
be 1000 km s _1 . We performed calculations with the 
mass-loss rates of M = 10~ 7 , 10~ 6 , 10~ 5 , 10~ 4 , and 
10 -3 M yr _1 . In the following, we especially focus on 
the two extreme cases, the models with M = 10 ~ 7 and 
10 -3 Mq yr _1 (hereafter they are referred to as the dense 
and dilute CSM models). 

The jet is injected from the inner boundary r = 10 9 
cm from t = to t = 60 s at a constant energy in- 
jection rate by using the same method as the previ- 
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Fig. 2. — Color-coded Lorentz factor and density distributions at 
t ~ 200 s for the model with M = 10 -3 M yr _1 (lower panel) 
and M = 10~ 7 Mq yr _1 (upper panel). 



ous works (e.g., IZhang et all 120031 : iMorsony et~aHl2QQ7t 
iMizuta et al.l l201lh . The parameters specifying the jet 
injection condition are as follows: the total energy E tot = 
3 x 10 52 erg, the energy injection rate E = 5 x 10 50 erg/s, 
the opening angle 6- } = 10°, the initial Lorentz factor 
Tq = 5, and the specific internal energy cq/c 2 = 20. 



3. RESULT 

3.1. Jet dynamics 

A lot of previous works on an ultra-relativistic jet ema- 
nating from the progenitor star have been carried out and 
unveiled the dynamical evolution of the jet, such as, the 
formation of the recollimation shock and t he realization 
of the well-known fireball solution, fe-g.. IZhang et ah! 
120031 : IMorsony et al.l [20071 : IMizuta et al.ll2Qllh . Our cal- 
culations successfully reproduce and confirm their find- 
ings. Some snapshots of the spatial distributions of the 
Lorentz factor and the density of the dilute CSM model 
are shown in Figure HJ The jet propagates in the in- 
terior of the progenitor star and then breaks out, and 
ejects stellar materials into the circumstellar space. As 
seen in the top right panel, the emergence of a hot ma- 
terial from the jet cavity follows the breakout of the 
collimated jet. The eject a rapidly expand to form a 
spherical cocoon as seen in the bottom left panel of 
Figure [TJ It is noteworthy that the cocoon expands at 
mildly relativist ic speeds. The appearance and the sub- 
sequent expansion of the cocoon have also been reported 
and investigated by several previous works (see, e.g., 
Alov et a l. 2000; Ramirez-Ru iz et al.ll200l IZhang et "all 
2003; Lazzati & Begelman 2005). 
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3.2. Effect of CSM interaction 
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Fig. 3. — Radial profiles along = 45° at t ~ 200 for the dense 
CSM model (solid line) and the dilute CSM model (dashed line). 
Each panel represents radial velocity normalized by the speed of 
light, the density, and the pressure from top to bottom. 



Results of the dense and dilute CSM models are com- 
pared in Figures [2] and [3] Figure [2] represents the spatial 
distribution of the Lorentz factor (left) and the pressure 
(right) at t ~ 200 s for the dense CSM model (lower 
panel) and the dilute CSM model (upper panel). 

Near the jet axis (6 < 10 — 20°), no difference between 
the two models is recognized. On the other hand, in 
the region with large inclination angles (6 > 20°), we 
can see differences between the models. Denser CSM 
reduces the size of the cocoon in comparison with dilute 
CSM. In addition, the pressure distribution shows shell- 
like structure in the dense CSM. 

This difference can also be seen in the radial profiles 
of some physical variables of the cocoon, as illustrated 
in Figure [3l In both cases, the expansion velocities are 
mildly relativistic as seen in the top panel. From the 
bottom panel showing the pressure profiles, one can see 
that the reverse shock forms as a result of the cocoon- 
dense CSM interaction. On the other hand, in dilute 
CSM, the rarefaction wave propagates toward the center 
in the cocoon. 

This is due to the aspherical distribution of the energy 
deposited into the ejecta. Near the jet axis, the energy 
carried by the jet is too enormous for the CSM to affect 
the jet propagation. On the other hand, the energy de- 
posited into the cocoon component is much smaller than 
that of the jet. The kinetic energy and mass of the co- 
coon component, which are now defined as those confined 
in the region outside the star and > 10°, can be ob- 
tained from results of the simulation. They are found to 
be 3 x 10 50 erg and 2 x 10~ 3 M at t = 10 s and 10 51 erg 
and 2 x 10 ~ 2 M at t = 20 s. These values are almost in- 
dependent of the mass-loss rate, because the dissipation 
of the kinetic energy of the jet to form the cocoon takes 
place in the star. The kinetic energy and mass of the 
cocoon increase due to the continuous energy and mass 



4 



SUZUKI&SHIGEYAMA 



injection by the jet. In our calculations, the injection of 
the jet is terminated at t = 60 s , which means that the 
injection of the mass and kinetic energy into the cocoon 
lasts even after the cocoon begins to expand. The kinetic 
energy of the cocoon is up to a few % of the total injected 
energy and thus has a potential for producing thermal X- 
ray photons with the observed luminosity ~ iq 45-48 erg 
s _1 for several hundreds seconds. On the other hand, the 
mass of the ultra-relativistic jet component, which is de- 
fined as the material with the Lorentz factor larger than 
100, is 3 x 10 -6 M , while a substantial fraction (~ 10 51 
erg) of the injected energy is carried by this component. 

As the radial profiles along 6 = 45° at t = 200 s in 
Figure [3] shows, a reverse shock is formed in the dense 
CSM model. The other model with the mass-loss rates 
10 -4 Mq yr _1 also form a reverse shock. The energy 
of the matter ejected immediately after the breakout of 
the jet from the surface results from the dissipation of a 
part of the kinetic energy of the jet for the initial several 
seconds. Denoting the dissipated kinetic energy by E^ 
and the fraction of the internal energy to the total by e, 
the pressure of the cocoon scales as 



47r(v ex pt) 



3 ' 



(2) 



where we have assumed that the cocoon is spherically 
expanding at the velocity v exp . On the other hand, the 
ram pressure of the CSM behind the forward shock is 
given by, 

MTL, 

(3) 
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/ c 2 ) 1//2 . The reverse shock forms 



when the pressure P c of the cocoon becomes comparable 
to the ram pressure p w T 2 c 2 of the shocked CSM. The 
balance between the pressure of the cocoon and the ram 
pressure yields the following expression for the time of 
the reverse shock formation, 



M^ex P r e 2 xp c 2 ~ U-0 x lO- 4 ^ Uo 51 erg 



(- 
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where we have derived the final expression by assuming 
^ex P = 0.9c. The value of the fraction e is found from 
the result of hydrodynamical calculations. This rough 
estimation is consistent with the fact that a reverse shock 
is observed for models with M = 10 -3 and 10 -4 M yr _1 
at t = 200 s and no reverse shock for models with lower 
mass-loss rates. 

3.3. Photospheric emission 

In the following, we investigate whether thermal X- 
ray emission from GRBs can probe their circumstellar 
environments. We derive the expected light curve and 
the spectra of thermal X-ray emission from our models 
by calculating th e photospheri c emis sion. 

According to iStarling et al.l (|2012[ ), thermal emission 
with the isotropic luminosity of the order of 10 4T erg 
s _1 and the photon temperature ~ 0.1 — 0.9 keV is ob- 
served. Illuminated by the radiation, heavy atoms, such 
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Fig. 4. — Light curves of the photospheric emission calculated 
for the models with M = 10 ~ 3 M yr _1 (solid line) and M = 
10 ~ 7 M© yr _1 (dashed line). 



as oxygen and carbon, are rapidly photo-ionized. The re- 
combination time scale much longer than the ionization 
time scale keeps those ions fully ionized and the domi- 
nant opacity source becomes electron scattering. Thus, 
we calculated the Thomson photosphere from a distant 
observer along the axis of the jet (6 = 0). In deriving 
light curves and spectra, we have assumed that the mat- 
ter and radiation are strongly coupled and the internal 
energy density on the photosphere is dominated by that 
of radiation. 

At first, we briefly consider the properties of the emis- 
sion. Since the ejecta move at mildly relativistic veloc- 
ities with the Lorentz factor of a few, the relativistic 
beaming effect strengthen the emission, especially, in the 
early phase. From the top and bottom panels of Figure 
[3j the radiation temperature of the shocked region for 
the dense CSM can be estimated to be, 



rfeBTp h = r/c B 



3p 



1/4 



0.1-0.2 keV, 



(5) 



in the observer frame. This is consistent with observed 
values. 

The resultant light curves and time-integrated vF v 
spectra of the photospheric emission for both models are 
presented in Figures 3] and At first, for both mod- 
els, the photospheric emission is very bright for the first 
~ 200 sec. This is because the cocoon is hot immedi- 
ately after the emergence from the stellar or wind photo- 
sphere. After the early phase, the radial velocity at the 
photosphere decreases as the photosphere moves inward 
and the cocoon gradually cools. This corresponds to the 
decrease of the luminosity. The cocoon cools in differ- 
ent ways for the dense and the dilute CSM. As seen in 
Figure [3l the cocoon component is shocked in the dense 
CSM model. The shock converts the kinetic energy of the 
cocoon into the thermal energy and keeps the shocked re- 
gion hot. As a result, the photospheric emission remains 
luminous even at t ~ 1000 sec. In the dilute CSM model, 
on the other hand, the cocoon adiabatically cools. 

Figure [5] shows the time-integrated vF v spectra for the 
dense (upper panel) and dilute (lower panel) CSM mod- 
els. In each panel, vF v spectra integrated over t = 0-200 
s (solid line), t = 200-1500 s (dashed line), and t = 0-1500 



THERMAL X-RAY EMISSION FROM LONG GRBS 



5 



Q. 



CD 



10 

10 
10 
10 
10 

10" 1 
10 
10 
10 
10 
10 

10" 





































































»■* 








































0-20 

— ?nnc 


Os 

-1500s 
00s 


























0-15 




















0.1 1 

photon energy [keV] 



Fig. 5. — Time-integrated vF v spectra of the photospheric emis- 
sion calculated for the models with M = 10 -3 Mq yr _1 (upper 
panel) and M = 10 -7 Mq yr _1 (lower panel). In each panel, uF u 
spectra integrated over t = 0-200 s (solid line), t = 200-1500 s 
(dashed line), and t = 0-1500 s (dotted line) are plotted. 



s (dotted line) are plotted. If we fit a blackbody spectrum 
with a single temperature to each of these spectra, the 
temperature is found to be k^T = 0.16, 0.077, and 0.13 
keV for the 0-200 s, 200-1500 s, and 0-1500 s spectra of 
the dense CSM model and k B T = 0.061, 0.038, and 0.051 
keV for the dilute CSM model. In high-energy part, how- 
ever, a deviation from the planck function is prominent. 
This shows that each spectrum is actually superposition 
of blackbody spectra with different temperatures. 

4. DISCUSSION AND CONCLUSIONS 

In this study, we performed hydrodynamical simula- 
tions of a jet emerging from a massive star surrounded 
by the CSM. Especially, we focus on the effect of the 
interaction between the eject a and the CSM. The CSM 
is assumed to be a steady wind with the wind velocity 
v w = 1000 km s _1 and the mass- loss rates ranging from 
M = 10~ 7 M Q yr _1 to 10 -3 M© yr _1 . We found that 
the dynamical behavior of the cocoon, which expands 
at sub-relativist ic speeds, is significantly affected by the 
ejecta-CSM interaction, while the collimated jet is so en- 
ergetic that the CSM can not decelerate it even for the 
dense CSM model. In the dense CSM, the cocoon is 
shocked and thus remains hot even at ^1000 s after the 
jet injection. On the other hand, in the dilute CSM, the 
cocoon cools adiabatically. Furthermore, calculating the 
photospheric emission from the cocoon, we found that 
the difference can be detected by observing their early 
thermal X-ray emission. 

Here we c ompa re our results with the calculation by 
iPe'er et al.l (|2006f ). They consider emission from a freely 
expanding spherical cocoon with the initial internal en- 
ergy of 3 x 10 51 , 10 52 , and 3 x 10 52 erg. On the other 



hand, the internal energy of the cocoon reproduced in 
the present calculations is less than these assumed val- 
ues. As a result, for the dilute CSM model, where a freely 
expanding cocoon is realized, the photon temperature of 
the photospheric emission is lower than 0.1 keV. For the 
dense CSM model, a fraction of the kinetic energy of the 
cocoon can be converted into the internal one by the re- 
verse shock, which leads to a higher photon temperature 

and brighter photosp heric emissi on. This effect is not 

taken into account by iPe'er et all (j2006l ) . 

From GRB 060218/SN 2006aj, a bright thermal X-ray 
emission with the temperature k B T ~ 0.1-0.2 keV is ob- 
served even in a few t housand seconds after the trigger. 
iCampana et al.l (j2006[ ) ascribed the long-lived thermal 
X-ray emission to a high mass-loss rate of M ~ 3 x 10 -4 
Mq yr _1 . Our results show that the shocked cocoon is 
realized at this mass-loss rate. Interestingly, some fun- 
damental features of the emission, such as the photon 
temperature and the long duration, are also reproduced 
by the photospheric emission from the shocked cocoon . 
Therefore, this event might occur in a dense circumstellar 
environment. However, we cannot reproduce the detailed 
temporal evolution of the observed photon temperature 
and the light curve by the present calculations in which 
the steady wind model is assumed. We attribute this 
discrepancy to the inhomogeneities of the CSM profile 
as described below. 

We should note that there are some uncertainties in 
this study. In particular, the properties of the photo- 
spheric emission calculated in this study are expected to 
strongly depend on the spatial distribution of the CSM. 
Although we adopted a steady wind model in this study, 
it may not be realized in actual circumstellar environ- 
ments of progenitor systems of long GRBs. Of course, 
the slope of the density profile of the wind depends on 
its mass-loss history prior to the gravitational collapse. If 
the progenitor had been rapidly rotating, the wind could 
have angular dependence. Furthermore, some authors 
point out that long GRB progenitors have evo lved in bi- 
nary systems (e.g.. iPodsiadlowski et al.1l2004f ). in which 
the structure of the material surrounding the system is 
expected to be much more complicated than single star 
progenitor cases. We regard the influence of the spatial 
distribution of the CSM on the properties of early ther- 
mal X-ray emission as one of future works. 
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